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THE PWN TORUS OF PSR J0538 + 2817 AND THE ORIGIN OF PULSAR VELOCITIES 
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ABSTRACT 

We find evidence for a faint wind nebula surrounding PSR J0538 + 2817 in CXO -ACIS imaging. This 
object is particularly interesting, as the pulsar spindown age is largest for any such X-ray PWN. If 
interpreted as an equatorial torus, the PWN supports the claimed association with the S147 supernova 
remnant and implies good alignment between the pulsar spin and space velocity. Comparison of the 
SNR age, X-ray cooling age and characteristic age suggests a birth spin period of > 130ms. In turn, if 
we accept as causal the alignment of the linear and angular momenta, this places strong constraints on 
the origin of the 'kick' at the neutron star birth. 

Subject headings: stars: neutron - pulsars - stars: rotation 
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1. INTRODUCTION 

Hubble Space Telescope and, especially, Chandra X- 
ray Observatory (CXO ) imaging have discovered toroidal 
structure in the synchrotron pulsar wind nebulae (PWN) 
of relativistic particles and magnetic fields around young 
neutron stars (e.g. Weisskopf et al. 2000, Pavlov et al. 
2001, Gotthelf 2001). Of these, the Crab and Vela pulsar 
tori are the most striking, with clear symmetry (pulsar 
spin) axes. Moreover, for these two objects HST proper 
motion studies show alignment between the position an- 
gles of the projected f2 and v. The chance probability 
of two such alignments is modest, but if taken as causal 
this gives strong constraints on nature of the 'kick' at pul- 
sar birth which typically produces a large space velocity 
(Spruit & Phinney 1998, Lai et al. 2001). Even a few fur- 
ther examples could improve the statistical significance of 
the Vt-v correlation and further constrain the kick physics 
(§4). The challenge is to identify sources for which both 
PWN tori and 2-D proper motions are measurable; after 
the Crab and Vela, many bright young PSR/PWN sources 
are quite distant and direct proper motion measurement 
is challenging. We report here on CXO imaging of one 
particularly interesting example. 

PSRJ0538+2817 , discovered by Anderson et al. (1996) 
is a 143ms pulsar with characteristic age t c — 6 x 10 5 y. 
Located within the old - 10 5 y SNR S147 (G180.0-1.7) 
with a DM-estimated distance of 1.2kpc (Cordes & Lazio 
2002) comparable to that of the SNR, it is believed to be a 
likely association despite the age discrepancy. Previously 
reported X-ray observations are limited to ROSAT all sky 
survey (Sun et al 1996) and HRI imaging. Inspection of 
the archival HRI data showed a source slightly extended 
beyond the nominal PSF, although aspect errors may have 
contributed. We have obtained ACIS -S imaging to resolve 
this possible PWN. 

2. CXO OBSERVATIONS 

We observed PSR J0538 + 2817 with the CXO ACIS 
on February 7-8, 2002, with the pulsar near the S3 chip 
aim-point, obtaining 19.5ks of TE exposure. All 6 ACIS 
chips were active and several field sources were detected off 
axis. None showed markedly unusual behavior; we there- 



fore concentrate on the S3 pulsar data. Calibration, pro- 
cessing and spectral fitting were performed with CIAO 2.3, 
including the latest calibrations and CTE corrections. 

2.1. Spatial Decomposition 

The pulsar core produces 0.13 cnts/s and, with the stan- 
dard 3.2s ACIS frame, suffers 20% pileup. In Figure 1 
we compare the 0.5-5 keV source radial profile with that 
of the HRMA+ACIS model PSF (generated for the ob- 
served source count spectrum and corrected for pileup). 
The point source normalization has been obtained from a 
PSF fit to the central 3". The core is resolved, with ~ 10% 
of the counts in extended emission. If the PSF is convolved 
with a a = 1" Gaussian, the counts are well matched out 
to 5", but a 5<7 excess remains above background from 
5 — 8". We thus believe that the extended emission cannot 
be attributed to aspect errors or poor PSF modeling. 
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Fig. 1. — Source and model radial count distributions. 

This emission does not follow the symmetric PSF wings 
(Figure 2). We start by fitting elliptical isophotes, with 
semi-major axes from 1-10", using the STSDAS routine 
ellipse. For small 1 — 2.5" semi-major axes, we find an 
cllipticity of ~ 0.1 — 0.2 at position angle ~ 130° (mea- 
sured N-E). From 4 — 9" the PA changes to a weighted 
average value of 61.0±5.5° with an ellipticity of 0.25-0.35. 
In Figure 2, the core has a tail to the SE defining the inner 
ellipse, while the scattered photons along the NE-SW axis 
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represent the larger ellipse. Some trailing along the read- 
out direction is visible in the detector coordinate image (at 
PA= 91°), but cannot be seen in Figure 2. The measured 
trail flux contributes less than 1 count in the 10" region of 
our extended emission, and does not lie along the axis of 
extension. 

In analogy to other pulsar/PWN systems we can ex- 
pect a toroidal wind shock surrounding the pulsar, pos- 
sibly with a polar jet. The 4 — 9" ellipse has a scale 
comparable to the expected wind shock radius for a con- 
finement pressure appropriate to the interior of S147 (see 
§3). We therefore adopt this physical model, computing 
the expected photons from the point source PSF wings 
and fitting the excess photons outside of 3" to an inclined 
torus. A uniform background for the residual counts is 
also allowed. The torus model is characterized by a radius 
Tt, a Gaussian blur torus thickness 5, the position angle 
^ of the torus axis and its inclination ( from the observer 
line of sight. We assume that the post-shock flow remains 
relativistic, so that the emission is Doppler boosted in the 
forward direction. For a rest power law emission spectrum 
dN' /dE^ oc E~ a , we expect an observed emissivity 

dN^/dE^ oc (1 - n • P)- {1+a) E- a 

where a is the photon spectral index, (5 = v/c is the bulk 
velocity of the post shock flow and n is the unit normal 
toward the observer (e.g. Polling et al. 1987). For typical 
PWNc, a w 1.5 and the emission brightens on the side 
facing the observer. 

Of the 2875 counts in the central 32" x 32" region, 
the best-fit point source provides 2430, the uniform back- 
ground contributes 80 counts and the 'torus' provides 50 
counts. The remaining counts represent unresolved fine 
structure in the nebula core. The best-fit torus radius 
is r T = 7.2 ± 0.7", the 'blur' scale is <5 = 1.9 ± 0.5". 
The torus (3 is 0.53 ± 0.04, similar to values fit for other 
PWN. The torus symmetry axis is located at position 
angle * = 154.0 ± 5.5° (measured N through E) and 
( = 100 ± 6° (into the plane of the sky), where the confi- 
dence regions have been scaled to that of ty. 
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Fig. 2. — Left: ACIS 0.5-5 keV image with smoothed low level con- 
tours. Right: Best-fit Point Source + torus model. Inset: Best-fit 
torus model at half scale, with <5 reduced to 0.5" to show geometry. 

In Figure 2 we show the photon data and the predicted 
image for the best-fit point source + torus model. The 
position angle W is robust to cuts in the fit data set and 
forced variations in tt and 5 and agrees well with the value 
determined from ellipse isophote fits; we adopt this value 
as determining the symmetry axis. £, primarily set by the 
ellipse axis ratio, is also surprisingly well determined. The 
brightening of the diffuse emission SE of the core indicates 



the near side of the torus, enhanced by Doppler beaming 
(inset). However, the core asymmetry also suggests a 'po- 
lar jet'. Certainly deeper exposures would be quite valu- 
able in confirming the reality of these faint structures and 
measuring their precise parameters. 

2.2. Spectral Analysis 

We extract the point source from a 2" aperture, correct- 
ing for pileup and measuring a nearby background from 
the S3 chip. For the torus we take a <~ 10" x 18" ellipse 
with a 4" radius exclusion core. The PSF wings of the 
central source contribute significant counts in this aper- 
ture, so we include a scaled version of the best-fit point 
source spectrum as a fixed background and fit the excess 
to characterize the torus emission. 
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Fig. 3. — Point source spectrum with a pileup-corrected magnetic 
H model atmosphere spectrum and residuals. 

The point source spectrum is adequately modeled by an 
absorbed blackbody, although systematic residuals persist 
at low energies. Table 1 gives the fit parameters with 
95% errors, the unabsorbed bolometric flux and the effec- 
tive radius of a spherical emittor at d = 1.2kpc. The 
second fit is from the magnetic H atmosphere models 
of Pavlov, et al. (1995), which produce slightly bet- 
ter x 2 and a larger effective radius. We note that all 
of these fits are somewhat suspect, since the spatial de- 
composition suggests <~ 10% diffuse emission even for the 
core. Composite (BB+BB or BB+PL) models of course 
improve the fit statistics, but are not demanded by the 
present data. These fits give Njj somewhat larger than 
the ~ 1.2x 10 21 cm~ 2 inferred from the pulsar DM with a 
canonical n e /n# ~ 10. 

A comparison of the energies of the 51 torus photons 
outside of 4" with the 240 photons from 1 — 2" (which 
avoid pile-up in the point source core) gives a KS statistic 
probability of 0.06 that they are drawn from the same dis- 
tribution. The torus region photons are thus significantly 
harder than the point source (despite the 30% contamina- 
tion by the PSF wings); this further supports the indepen- 
dant nature of this component. We assume a PL model 
for this torus emission. Fixing Nh as above and assum- 
ing a typical PWN T = 1.5, we find an unabsorbed flux 
/x(0.5 - 5 keV) w 1.6 ± 0.4 x 10- 14 erg/cm 2 /s for the full 
torus. This spindown efficiency r\ w 6 x 10~ 5 d 2 2 is about 
10 x lower than that of most r <~ 10 4 y pulsars. 
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3. THE SNR ASSOCIATION & KICK-SPIN CORRELATION 

S147 is a late blast wave phase SNR, quite circu- 
lar in both radio (Furst & Reich 1986) and Ha (van 
den Bergh 1978) images, with a ~ 166' diameter. 
The estimated 0.8-1.6 kpc distance agrees very well 
with the pulsar's 1.2 kpc DM-derived distance. Sofue, 
Furst & Hirth (1980) find the geometrical center at 
05:40, +27:44 (±2.5') (J2000). At the nominal distance 
the 7 = 5/3 Sedov solution gives a blast wave age of 

t = 1.4 x lO 5 ^ 2 (ra/£5i) 1/2 y and an expansion velocity 
v cxp = 82d~ 3 2 /2 (E 51 /n) 1 / 2 km/s for a local ISM density 
n cm~ 3 and an explosion energy 10 51 i?5ierg. The op- 
tical expansion velocity 80km/s (Kirshncr & Arnold 
1979) agrees well, supporting the distance estimate. 




Fig. 4.— llcm image of S147 (after Furst & Reich 1986), with 
PSR J0538 + 2817 's proper motion vector (3° field, N up/E left). 

If we assume magnetic dipole spindown (n — 3) with a 
constant field 10 12 £?i2G, the spindown luminosity at char- 
acteristic age 10 5 T5y is 

E= 1.0 x 10 36 (Bi 2 r 5 )- 2 erg/s, 
leading to a scale for the wind termination 

r T « {E/AKcP^t) 1 ' 2 . 
This termination will be azimuthally symmetric (e.g. 
toroidal) whenever the external static pressure P ex t > 
Pram = 6 x 10~ 10 m; 2 g/cm/s 2 , i.e. subsonic pulsar mo- 
tion at a speed lOOi^km/s. In the general ISM, P ext 
is typically a fewx 10~ 13 g/cm/s 2 ; all but the slowest 
young pulsars should display non-toroidal bow shocks. 
However, in the interior of a Sedov-phase SNR we have 
P ext ~ 6 x 10~ 11 (n^| 1 /^) 1 / 5 g/cm/s 2 where the ambi- 
ent ISM density is n and we have distinguished between 
the true (=SNR) pulsar age t and the characteristic age 
t c . Since the density in the hot interior is quite low, 
eg. Log(rii) = Log[no(2r/Vsod — 2)] (Mansfield & Salpeter 
1974), we expect pulsars in the interior of SNR to be in 
the toroidal shock phase for t 5 < 3. The difficulty is that 
at typical speeds of v-j 5, a pulsar outruns the blast 
wave by t 5 ss 0.3(E 51 /n ) 1 / 3 (v 7 /5)^ 5 ^ 3 . In fact, as soon 
as it passes the contact discontinuity, the large density and 
ram pressure jump ensure a bow shock morphology. PSR 
1951+32 in CTB 80 is likely in this phase (Shull et al. 
1988). 

According to these estimates, a torus for PSR J0538 + 
2817 is not unexpected. At t$ > 1 we would normally find 



it in the ambient ISM, but its low 1.5/ts keeps it in 
its host SNR. With a low Bi 2 = 0.74, the pulsar husbands 
its spin-down luminosity giving a torus of angular size 

Otor = rr/d w ^"(Biadkpc-nsng-^tB °- 6 )- 1 

or about 3"t° 5 6 n °- 3 d^ 1 2 for PSR J0538 + 2817 , among 
the largest angular sizes expected for known pulsars (after 
Crab and Vela). Since the observed inner edge of the torus 
is at ~ 5", these estimates agree if t 5 ~ 2 or ii ~ 0.2. 

Given the large offset of PSRJ0538 + 2817 from the rem- 
nant center, a good test of the Q-v correlation can be made. 
The offset vector implies a proper motion (23.6 ± 1.5) /t^ 
mas/y at position angle —32.3+ 3.6°. The pulsar trans- 
verse velocity is 134di. 2 /i5 km/s. Together with the torus 
geometry estimates of §2.2, this offset gives a kick-spin po- 
sition angle difference of Qq,- v = 6.3 ± 6.6°. If v is along 
the torus axis, the 3D velocity is 137di.2/i5km/s, rather 
low for an isolated pulsar. 

The Crab pulsar optical proper motion (Caraveo & 
Mignani 1999) and inner nebula axis vectors give 9q~ v = 
10 ± 10°. For Vela, the symmetry axis of the X-ray torus 
in Pavlov, et al. (2001) is aligned within 2 + 3° (Bailes 
et al. 1989) or 10+3° (De Luca, Mignani & Caraveo 2000) 
using the radio or optical v, respectively. The joint prob- 
ability of obtaining these (2-D) projections by chance is 
~ 0.3 — 1%, making the association intriguing, but not 
definitive. If one adopt PSR J0538 + 2817 as a third align- 
ment, the chance probability drops to <~ 2 x 10~ 4 . 

4. DISCUSSION: AGE AND KICK CONSTRAINTS 

One of the key discriminants of the origin of this corre- 
lation is the initial pulsar spin period Pq (Lai, Chernoff & 
Cordes 2001). For a constant braking index n one finds 

t = 2T c [l-(P /P)("- 1 >]/(n-l). 

Since our PWN position angle supports the connection 
with S147, we must reconcile the r c = 10 5 ' 8 y (n = 3) char- 
acteristic age with the maximum plausible age ~ 10 5 y 
of the SNR. One possibility is very rapid B decay, giv- 
ing an effective n > 10. More likely, as for several other 
PSR-SNR associations there is a large initial spin period. 
The observed thermal luminosity supports the young age 
and large P . This L w 2 x 10 33 <f 2 2 erg/s, interpreted 
as full surface cooling agrees well with the flux in stan- 
dard cooling scenarios near 10 5 y. At Pq rs 130ms (is w 1) 
PSRJ0538+2817 has one of the slowest initial spins known. 

Physical kick-spin connections generally rely on rota- 
tional averaging (Spruit & Phinney 1998), so this large Pq 
is problematic for many proposed mechanisms. If the kicks 
are isotropically distributed, the net moment at 9 w 60° 
must be averaged for Tkick S> 0.38-PoA an $o-t> to achieve 
a final spin-kick angle 6q- v . Indeed, since 'Natal' kicks 
may be applied before the neutron star has reached its 
final ~ 10 6 cm radius, conservation of angular momen- 
tum for kicks at r$ > 1 increases the effective initial 
period by r 2 . Thus for PSR J0538 + 2817 , we require 
Tkick > 0.5 s(l — t5/6) 1 / 2 r| for the observed 0q- v . 

Certainly the E-M rocket (Harrison- Tademaru) effect 
is hopeless, as it requires Po < 3ms for lOOkm/s veloc- 
ities. Also hydromagnetic models imparting momentum 
over ~ 0.1s are inadequate, especially as the momentum 
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deposition is at the re > 10 shock break-out. Even mech- 
anisms invoking asymmetric v emission are severely con- 
strained, but appear just feasible considering the ~ 3s v 
diffusion timescale at r 6 <~ 1 — 2. However inducing the 
large required v asymmetry is challenging. For magnetic- 
field induced asymmetry, even PSR J0538 + 2817 's mod- 
est velocity requires an ordered neutrinosphere field of 
<~ 10 15 5 G (Arras & Lai 1999), otherwise unmotivated 
for this low B pulsar. 

As an alternative one may imagine schemes where the 
spin is driven to the kick axis. One possibility is post- 
collapse accretion of a remnant disk (Blandford, private 
comm.). Imagine a disk of ~ AM accreting at a super- 
Eddington rate onto a 1.5Mi^Mq neutron star with mo- 
ment of inertia 10 45 /45gcm 2 . The disk is approximately 
Keplerian down to some inner radius r in , at which point 
we assume roughly half of the mass accretes while the 
remainder is ejected in polar jets with velocity v esc w 
2 1//2 VKep('"in) an d an asymmetry \. The large Pq in- 
ferred for PSR J0538 + 2817 is a serious constraint on 
this scenario. If the disk extends to the star, we might 
expect jet velocities of ~ c/4 as in SS443, but only 
AM w 2 x 10 _3 /4 5 (100ms/Po) M may accrete before 
one spins up to P < 0.1s. With this mass, jet recoil gives 
a kick v K w 75x/M 1 . 5 km/s for PSRJ0538 + 2817 's param- 
eters, so that even a large jet asymmetry does not produce 
the desired momentum. 

Perhaps more realistically, standard equilibrium spin- 
up arguments set r in = (B^/My/2 G M*) 2 / 7 w 1.8 x 



10 8 ' M_^"' cm. The Keplerian period at this Alfven 

radius is thus comparable to P for M_§ ~ 130, i.e. 
~ 100 x the Eddington accretion rate. If as expected the 
star is born near a ~ 1 P^ms break up period, the mini- 
mum mass needed to spin down the star to P is AM w 
0.04M Q / 4 5Pr lM iT5 /3 (- p o/100ms)- 1 / 3 , which gives a kick 

velocity v k w 800x/45M 1 " 5 1/3 (P /100ms)- 2 / 3 /P km/s. 
PSR J0538 + 2817 is relatively slow with v w 140 km/s 
and Pq « 130ms, so we require an jet asymmetry \ w 

1/3 

0.2PiM^/I 45 . Since the velocities are set at r in we see 
that slow spin goes with slow kick in this model. Sub- 
stantial jet asymmetries are required, but remain to be 
explained. Of course, additional mass accreted while at 
equilibrium reduces the required x- 

With its large P and small n _ v , PSR J0538 + 
2817 provides some of the strongest constraints on spin- 
kick physics. These important physics conclusions are 
presently based on rather limited data; clearly better im- 
ages confirming the torus and refining W are needed. Also, 
a radio interferometric proper motion can improve the ve- 
locity vector and get an independent kinematic measure 
of £5 . Such data should provide a new window on neutron 
star core collapse and its immediate aftermath. 
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Table 1 
Point Source Spectral Fits 

Model kT(keV) N H (HFcm" 2 ) / x (10- 12 erg/cm 2 /s) P e //(km) X 2 /DOF 



BB 0.159±0.0017 3.1±0.2 3.3 2.6 1.03 

magH 0.054±0.0007 4.7±0.2 13.7 13 0.94 



